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Abstract A new expression for the characteristic function of log-spot in Heston model is presented. 
This expression more clearly exhibits its properties as an analytic characteristic function and allows 
us to compute the exact domain of the moment generating function. This result is then applied to the 
volatility smile at extreme strikes and to the control of the moments of spot. We also give a factor- 
ization of the moment generating function as product of Bessel type factors, and an approximating 
sequence to the law of log-spot is deduced. 
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, , 1 Introduction 

Ph The first surprising fact about the Heston stochastic volatility model (Heston [1 1]) is that the charac- 
teristic function of log-spot is computable and has a nice expression in terms of elementary functions; 
its deduction had enormous merit. The second thing, and still more fascinating, is that such character- 

^ istic function is analytic, that means (see Lukacs [16, chapter 7] for an equivalent definition and the 
main properties of analytic characteristic functions) there is a function ^(^) of the complex variable 
^ z, analytic in a neighborhood of 0, such that 

for u (real) in a neighborhood of 0, where Xt is log-spot at time t. 
(3 The fact that a characteristic function is analytic has important consequences, and we here exploit 

some of them. However, the standard form of writing the Heston characteristic function hides the 
analycity property and difficults its utilization. So we first give a new expression of the characteris- 
. ^ tic function, but more importantly, we use that the analycity property is equivalent that the random 
^ variable has (real) moment generating function in a neighborhood of the origin: there is £ > such 
S that 

M(u) = E[e"^*] < oo, -e <u<e. 
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In that case, the function M{u) is a real analytic function in (— e, e) and the main properties of the 
characteristic function can be studied through M{u), which is a real function, and in general much 
simpler to analyze. 

As a consequence of the study of the moment generating function we obtain the domain of that 
function and we give a numerical simple procedure to compute the poles of the characteristic function 
in its strip of convergence. This has several practical consequences and we apply it to the computation 
of the smile wings parameters in the formula given by Lee [14]. We also apply these results to the 
assessment of the moments of a particular model; such study is complementary of the one of Andersen 
and Piterbarg [2]. 

The second contribution of this paper is that we factorize the moment generating function of log- 
spot. This allows us to identify a Bessel type densities as the building blocks of the log-spot. This 
result has interesting applications. For example, it allows us to construct a sequence of random vari- 
ables converging in law to log-spot. On the other hand, for a certain combinations of the parameters, 
log-spot Xt is a sum of independent non-centered random variables, and we can identify Xt as 
a member of the non-homogeneous second Wiener chaos (see Janson [13, chapter 6]); this agrees 
with the intuition that for some parameters, the stochastic volatitlity, given by a CIR model (Cox et 
al. [3]), has the law of the sum of the squares of a finite number of independent Omstein-Uhlenbeck 
processes, that are in the second Wiener chaos, and such property is transferred to Xt. 

The paper is organized as follows. First we deduce the moment generating function and the char- 
acteristic function of log-spot; as far as we know, these expressions are new. Both can be obtained 
manipulating the expressions obtained by Dufresne [6] or the standard expressions of the characteris- 
tic functions (see for example, Gatheral [9] or Albrecher et al. [1]), however we prefer to give a new 
deduction from scratch since our procedure is quite general and can be applied to other problems. In 
Section 3, we obtain the domain of the moment generating function. In Section 4, we give some appli- 
cations. Finally, in Section 5, using techniques of complex analysis, the moment generating function 
of log-spot is factorized. The analysis of the factors allows us to identify them as moment generating 
functions of Bessel type random variables, and to construct a sequence of random variables that con- 
verges in law to log spot. In the Appendix we review some facts on moment generating function of a 
random variables and we put technical details of the proofs. 



with initial conditions Sq = sq > Q and Vq = vq >Q, where a, 6 > and cgM — {0}are constants, 
and W and Z are two standard correlated Brownian motions, {Z, W)t = pt, for some p e [—1, 1]. 
The process is a Feller diffusion (Feller [7]) or, in the financial literature, a CIR model (Cox et al. 
[3]). The parameter a is called the the mean reversion factor, b is called the long term volatility and it 
is also written V^o, and c is called the vol-of-vol. Write 



2 The Heston model 



The Heston model [11] is defined by the system of stochastic differential equations 




(1) 



Xt = log St- pt. 
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By the Ito formula, 

dXt = -]^Vtdt+^tdZt, 
with initial condition Xq — xq — In sq. Thus, we will consider the system 

dXt = --Vtdt+ JVtdZt 1 
dVt = a{h - Vt) dt + c^/VtdW{t) J 

2.1 The moment generating function of Xt 

First, we check that Xt (indeed {Xf, Vt)) has moment generating function, and later we deduce its 
expression as the solution of a (real) PDE obtained from Ito formula. 

For every v e M, the random variable 6^"^*+^^* is positive, so we can compute its expectation 
but it can be infinite. Write 

Zt^pW{t) + y/i-pm'(t), 

where W is a Brownian motion independent of W, and use the habitual trick 

g{u,v,t) := E[e"^*+^^*] = E[E[e"^*+^^7iy'(s), < s < 

We obtain 

( uvqP upabt-. 
g{u,v,t) = exp [xqu —j 

E 



exp { (" + ^) V; + (y - 2 - V + ^) X ^' 



Note that the coefficient of Vg is the equation of a parabola in u through the origin, so if u is near 

zero, that coefficient should be also near zero. Since both Vt and Vs ds have moment generation 
function (see Dufresne [5]) it follows that the expectation is finite for {ua;) in a neighborhood of 
(0, 0). Fix T > 0. Applying the Ito formula to exp{MXj + vVt] and taking expectations, we have 

g{u, V, t) = e""°+'"'« + vab f g{u, v,s) ds + ( - + + + uvpc - va) f ^^^''^'"^^^ ds, 

Jo \ 2 2 2 ■'Jo tjv 

where we have used the property of the moment generating function 

dv 

Differentiating with respect to t, we get 

dg(u,v,s) dg(u,v,t) 

0^ P{U:V) — = abvg{u,v,t), 

where 

p{u, v) = -- + — + — + uvpc - va. 
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This equation has a unique solution that is 

(lb I (? 

( .\ [p{'^,<f>{u,v,t))\ . abtpu a?ht. 

g{u, v,t)^ r exp [uxo + (l){u, v, t)Vo \ ^ | , 

\ p{u,v) I ^ c & ' 

where 

P{u) — — pcuY + c^{u — v?)^ 
^{u,v) — —2 arctanh((c^v + cpit — a)/P(it))/P(it) 

^{u, v,t) = + ^ - ^P{u) tanh (P{u) {t + 7(m, v))/2^ . 

For = 0, we get the moment generating function of Xt, Mt{u); when there is no confusion we will 
suppress the subindex t and write M{u). After some tedious manipulations, M{u) can be written as 

M{u) =E expjwXt}] 

^{a-cpu)t/2 \ 2ab/c^ 



exp{xQu} 



cosh{P{u)t/2) + (a - cpu)smh{P{u)t/2)/P{u) 



{u-u^)smh{P{u)t/2)/P{u) \ 
' 1 cosh{P{u)t/2) + (a - cpu) smh{P{u)t/2)/P{u) J ' 

where 

P{u) = \/ {a — pcuY + c^{u — V?). 

Remark 2.1. Formula (3) coincides with the one that can be deduced from the joint Laplace-Mellin 
transformation of S't, 1^, Vg ds given by Dufresne [6, Theorem 12]. 

2.2 The characteristic function of Xt 

For z complex, consider the function 

^{a-cpz)t/2 \ "iab/c^ 



cosh(P(^)t/2) + (a - cpz) smh{P{z)t/2)/ P{z] 



r_ {z-z^)smh{P{z)t/2)/P{z) \ 

■ ^""^ \ cosh{P{z)t/2) + (a - cp;^) sinh(P(2)t/2)/P(z) J " ^ ^ 

Writep(2;) = (a— pc2;)^ + c^(z — z^) the second degree polynomial within P (2;). Sincep(O) = > 0, 
using standard techniques of complex analysis, we see that ^{z) is well defined and analytic in a 
neighborhood of 0. Obviously, = M{u) on a (real) neighborhood of 0. Then (see the Appendix, 
Proposition A.l, and Section 5), the characteristic function of Xt is Explicitly, for e M, 

/ gat/2 \ 2o6/c2 

(p(u) = E[e^"^*] = = explixowj — , , , , , ,,,,,,, 

^ ^ ^ ^ ^ ^'^ " ^ \cosh{dt/2) +^smh{dt/2)/dJ 

j {iu + u^)smh{dt/2)/d 1 

• exp I - vo +^smh{dt/2)/d)' ^ ^ 
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where 

d = d{u) = P{iu) = \/ {a — cpiuY + c^{iu + v?)^ 
^ = ^{''A — ^ ~ cpui. 
After some computations we arrive to the formula of Albrecher et al. [1] 

ip{u) = exp{ixou} exp I - (^(e - d)t - 2 log j | exp | - - (i) ^Zr^^t / ' 

where 

9 = 9W = 

i + d 

Of course, formula (5) looks more complex than the compact (6). However, when one recovers from 
the shock, one realizes that the former is easier to handle that the latter. 



2.3 Inversion of Heston model and some comments on the parameters 

We will say that a process S = {St, t G [O.T]} is a Heston type process, and we write S ~ 
'HVv{a, 6, c, p, So, vq, jj) if St verifies a system (2) with some stochastic volatility Vt. When there is 
no confusion with the underlying probability P we will omit it. 

Recent results of del Bano [4] show that if 5" ~ HVF{a, b, c, p, sq, vq, ji), a > cp, then 

5"^ ~ HQ{a - cp, ab/{a - cp), c, -p, \ Vq, -p), (7) 

where Q is the probability given by 




dF 



To prove that property it is needed to work with the whole process S. However, an easy verification 
can be done using the moment generating function (3). We will use this property in Section 3 to 
compute the values of E[exp{MXt}] for negative values of u. 
It can also be proved that 

HVwia, b, c, p. So, vq, p) ~ HV^{a, b, -c, -p, sq, vq, p). 

Again, it is needed to consider the process to prove this equality, but a check is deduced from the 
expression of the moment generating function (3). So, without loss of generality we will assume 
from now on that c > 0. 



3 The domain of the moment generating function 

In this section we deduce the domain of the moment generating function (3); this deduction is not 
direct since we obtained M{u) not by the computation of the expectation E[e"^*] but by an indirect 
way. So, we only know that the moment generating function coincides with the function given in the 
right hand side of (3) in a neighborhood of zero. However, as stated by Lemma 3.3 below, since the 
function is analytic, we are in safe land. In the first subsection, using this idea we do a first study of the 
moment generating function. Later, in the second subsection we work with the analytic continuation 
of the function in the right hand side of (3). 
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3.1 Preliminary study of the domain of the moment generating function 

The main ingredient of the moment generating function given in (3) is 

f{u) := cosh(P(«)t/2) + (a - cpu) ''""^ pff^ , (8) 

where 

P{u) — s/ {a — pcuY + c^{u — V?). 

Write 

p{u) := (a — pcuY + (?{u — u^). (9) 

When p 7^ ±1, p{u) represents a parabola with leading coefficient c^(p^ — 1) < 0, and for w = 0, we 
have > 0. So, it is an inverted parabola with real roots, U- <Q < u+, given by 



c — 2ap ± \/^a? + — 4acp 

^ MI^T^ ■ ^^'^ 

When p——l, then 

= + c(c + 2a)ii, 

is a straight Une with positive slope that intersects the horizontal axis at — —a? / (c(c + 2a)) < 0, 
and we write u+ — oo. 

Similarly, for p — 1, p{u) degenerates in the straight line 

p{u) = + c(c — 2a)u, 

and the slope can be negative, postitive or zero or negative, and 

• When 2a > c, then — — oo. 

• When 2a < c, then u+ = oo. 

• When 2a — c, then u- — — oo and u+ — oo. 

Thus, for every p E [—1, 1], we have that p{u) > on [u-,u+], hence the function f{u) is well 
defined and analytic in such (possible infinite) interval. Denote by D{Xt) the domain of M{u), 

D{Xt) = {« e R : M{u) = E[e"^'] < oo}. 

In principle (see Lemma 3.3 bellow) M{u) is defined in the subinterval of between the 

biggest negative zero and the smallest positive zero of f{u). Next proposition summarizes the study 
of such zeroes in that interval. 

Proposition 3.1. For every t > 0, and p e [—1, 1] there is the inclusion [u-, 1] C D{Xt). Moreover, 

1. When a > cp ( in particular, for every p < 0), then for all t > 0, the function f {u) has no zeroes in 

u+], and consequently, u+] C D{Xt). (Except for p = I and c = 2a). 

2. When a < cp, write to — 2/{cpu+ — a) >0 

(i) Ift < to, then [u-, u+] C D{Xt), and f{u) has no zeroes in [u-, 
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(ii) Ift > to, then f{u) has no zeroes in [n , 1], and has one and only one zero in (1, 
3. When p=landa = 2c, then D{Xt) = ( - oo, 1/(1 - e""*)). 

Remark 3.2. 

1. The cases p = ±1 are specially important. We stress that 

(i) For p = — 1, we are always in case 1 and = oo. 

(ii) For p = 1, when a > c, we are in case 1 (except if a = 2c). For a < c we are in case 2; 
however, if a < c < 2a, then u+ = oo, and we have to = 0, and thus we are in case (ii) for 
all t > 0. 

2. From the preceding proposition it follows that E[5't] < cxo for every a,b,c > and p G [—1,1]. 
Moreover, when p = 0, by construction, {St, t E [0,T]} is an exponential local martingale, that 
is a positive supermartingale, see Revuz and Yor [17, pages 148 and 149], and 

E[St] = M{l)=e^'>, 

so {St, t E [0, T] } is a true martingale. This was proved by Andersen and Piterbarg [2, Proposition 
2.5] using the Feller explosion criteria and Girsanov Theorem. 

3. As a continuation of the preceding point, we should remark that {St, t E [0, T]} is not always a 
square integrable martingale. This can cause some problems. See Subsection 4.2. 



3.2 Computation of the abscissae of convergence of the moment generating 
function 

When there is no root of f{u) in [1,m+] the domain D{Xt) is larger than (in particular, 

when a > pc). To carry out this study, we need more properties of the moment generating function. 
Consider an arbitrary random variable X, with moment generating function Mx{u), and domain 
Dx = {u eR : Mx (u) < cxo}. Remember that Dx is a interval of IR (finite or infinite, open or 
closed from one side or the other, that always include the origin, and it may be just {0}), and Mx 
is analytic in the interior of Dx- The left (respectively right) extreme of Dx is called the left (resp. 
right) abscissa of convergence, and plays a major role. The following property is well known, but 
since it is key in this paper, we stress it. We give the property for the right abscissa of convergence, 
and a similar statement is true for the left abscissa. 

Lemma 3.3. Let X a random variable such that there is a neighborhood of zero included in Dx- 
Assume that (r, s] C Dx, and that there is an analytic function h : (p, g) ^ M such that 

1. (r, s] C (p, q). 

2. Mxiu) = h{u),foru E (r, s]. 

Then Mx = h on (p, q). Moreover, if lim„^g h{u) = oo, then the right-abscissa of convergence of 
Mx is the point q. 
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Proof. 

Denote the interior of Dx by {a, 13). li (3 = oo, by analytic continuation, Mx = h on (r, g). 
Consider /3 < oo, so /3 is the finite right-abscissa of convergence. Then the function Mx has a 
singularity at P (see the Appendix, Proposition A.2) Thus, s ^ (3, because Mx{s) = h{s) and h is 
analytic in s. Hence, (3 > s. In the same way, (3 < qh contradictory, then (3 > q and by analytic 
continuation, Mx = hon (r, q). The second part of the Lemma is obvious. □ 

Now we apply the preceding lemma to the moment generating function of the log-spot, M{u). 
When the function f{u) given in (8) has no zeros in [1, u+\, since 

lim M{u) < oo and lim M(u) < oo, 

u\u- u^u+ 

by Lemma 3.3, the domain of M{u), D{Xt), is bigger than u+] Then, to assess more carefully 
that domain, consider the function cosh ^/x, for x > 0. Its Taylor expansion is 



cosh ^/x — ^2 



X" 



n=0 



(2n)!' 



The series on the right defines an entire function, say Li(x). However, when x < 0, that series 
coincides with the Taylor expansion of cos ^/—x. Hence, Li{x) is an entire function that, when 
written as the composition of elementary functions, has different expression according to whether 
X > or a; < 0, that is 

cosh y/x, if X >0, 



Li{x) 



cos \/—x, if a; < 0. 



In a similar way, the function (sinh y/x)/y/x, x > Q, can be analytically continuated to negative 
numbers as (sin a/^) / y/—x, x < 0. 

Denote hyu*_ > — oo the left abscissa of convergence of M{u) and by < oo the right abscissa. 

Put 

P{u) = V-P{u). 

For u e {u*_,U-) or M e {u+, u*_^), by Lemma 3.3 the moment generating function is 

^{a-cpu)t/2 \ 'iab/c^ 



M{u) — expjxoM} 



cos(P(M)t/2) + (a - cpu)s\n{P{u)t/2)/P{u) 



exp ^ - ^ (., -«')sm(P(„)t/2)/P(„.) 

cos{P{u)t/2) + (a - cpu) sm{P{u)t/2)/P{u) 



In that expression, the main part is the function 



f{u) := cos{P{u)t/2) + {a- cpu f^^S''^^''^\ (12) 

P{u) 



Heston characteristic function 



9 



Both / and / are defined in disjoint sets, and can be combined in a new function 

{f{u), Hue [u-,u+], 
(13) 
f{u), if u < or u > u+, 

that is analytic in R. See Figure 1 for a plot of that function. 



/ X X X u*/u 






X X X / • 





Figure 1. Plot of F{u). The points u*_ and u*_^_ are the abcissae of convergence of M{u). 

Then, to find the right abscissa of convergence, we need to find the zero of the function f{u) in 
[1, and if there is no zero, then to find the smallest zero u > u+ of f{u). Note also that the real 
zeroes of f{u) are the real solutions of the equation (see Figure 2) 

tan(P(ii)i/2) = (14) 

a — cpu 

(Except when a/cp — k7r/2 for some natural number k). 





1 i 
1 / 


1 

r 

/ 

\ u* 
\ — f 


\ 




w 


1 





















Figure 2. Solid line: Plot of tanh {P{u)t/2), for u e or tan [P{u)t/2) in the 

complementary. Dashed line: Plot of —P{u) /{a — pcu) or —P{u) /{a — pcu). 

For the left abscissa of convergence we need to look for the biggest solution m < m_ of f{u) = 0, 
or, equivalently, to work with equation (14). 
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In order to give a bound for the abscisses of convergence, denote by a±i the solutions of the 
equation 

47r2 

that for all t > are real and a_i < u- and a+i > m+. Also put /3±i the solutions of p{u) = — 27r^/t^. 
Consider the following cases: 

1. If a > pc, then u\_ E a+i). This can be deduced from the consideration that the image of the 
function ta,n{P(u)t/2) on that interval is M, and the properties of the function in the right hand 
side of (12). The only case not clear is when a/{cp) = tt/2, due to the fact that —P{u) /(a — cpu) 
has a vertical asymptote ?lX.u — Pi, this case is studied by direct inspection. 

2. If a < pc, then remember that / (I) > 1, and, on the other hand, 

/(/3+i) = (a-cp/3+i)2/7r<0, 
because f3^i > 1 and p > 0. So F{u) has at least one root in (1, P+i). 
Joining these comments with Proposition 3.1 we have 
Theorem 3.4. With the above notations, 

1. If a > pc, the right abscissa of convergence is the smallest zero u > of f{u) — 0, and 
u'^ e {u+, a+i). (Except for p — 1 and a — 2c). 

2. If a = pc, the right abscissa of convergence is u'^ — 1. 

3. If a < pc, let to = 2/ {cpu^ — a) > 0. 

(i) Ift < to, then u'^ is smallest zero u > u+ of f{u) — 0, and u'^ e {u+, 

(ii) Ift>to, then u*^_ is the zero of f{u) in {1, u+]. 

4. In every case, the left abscissa of convergence, u*_ is the biggest zero u < U- of f{u) = and 
ul_ e U-). 

Remark 3.5. 

1 . When a > pc the inversion formula (7) can be used to compute u*_ in terms of u*_^_ of the inverted 
model. Specifically, 

u*_{a, c, p, t) = — M^(a — cp, c, —p, t) + 1. 

2. This theorem gives a direct procedure to invert the formulas of Andersen and Pitebarg [2, Propo- 
sition 3.1]. 

4 Applications 

In this section we present some applications of the exact knowledge of the abcissae of convergence of 
the moment generating function of log-spot. 
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4.1 The smile at extreme strikes 

One of the motivations for the study of the abcissae of convergence of the moment generating function 
of log-spot in the Heston model in this work is the outstanding result of Roger Lee ([14]) where an 
explicit relation is found between the assymptotic behaviour of the volatility smile and these abcissae. 
This can be of interest in designing sensible smile interpolation and extrapolation schemes as has 
been shown by Gatheral [8]. In the case under consideration, that of the Heston dynamics, if u\_ and 
u*_ are the abcissae of convergence of log-spot, then according to Lee [14] the asymptotic behaviour 
of the Heston smile for expiry T as the strike K goes to infinity is o-{K) ^ l3R\n[K)/T where 
(5r e [0, 2] is defined by 

Likewise the behaviour as K approaches zero is (t{K) ^ y/—(3L^'n.(K)/T where (3l e [0, 2] is 
defined by 

2/3^+8 2 " 

As an example, consider market data for a one year equity smile p — —90%, a — 2, c — 80%, b — 
15%^ and t — 1. The calculations outlined above yield u*_^_ — 37.43 and u*_ — —3.21. By Roger Lee's 
formulas above this implies the behaviour of the smile wings is ruled by 

Pr = 0.01 and Pl = 0.13 

where, as expected, the smile at low strikes has a larger coefficient. These numbers can be of use to 
design an extrapolation scheme for the Heston smile in extreme strikes where the numerical integra- 
tion breaks down. 



4.2 The importance of second order moment of spot in stochastic volatility 
models 

The moments of spot and the abcissae of convergence of log-spot are related because 

< = sup{u e K : E[S^] < oo} and u*_ = mi{u e M : E[S^] < oo}. (15) 

So, if u*^_ < 2, then the second order moment of spot is infinite. This can cause problems in pricing 
certain standard European derivatives, as the next couple of examples show. 

1. Pricing an FX performance note. Consider a performance note that pays Euros on the perfor- 
mance of the EUR/USD exchange rate. This is a contract that pays the following amount in Euros at 
expiry T 

Fayoff (EURUSDt) = Notionabc,^ [^^^^ - 1^ EUR 

The notation should be self explanatory: EU RU SDt is the EUR/USD exchange rate at time t and 
t = is today. By the standard martingale methods the price of a derivative product is simply its 
discounted expectation, this might (and does) lead naive market participants to price such a transaction 
as 

Present Value = L'/I^^^Notionabc;,^ { eURUSDo ~ 
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where we have used the fact that the risk neutral expectation of spot is the forward Ft, and write 
Df^^^ for the relevant discount factor. Unfortunately this approach is flawed since EURAJSD is the 
price in US Dollars of one Euro and in the equations above we are basing a Euro payment on this 
quantity. Its correct price can be derived by noting that its payout in Dollars is 

FayoffiEURUSDr) = NotionalEUR ( ^.^fj^ - 1 ) ■ EURUSDt USD 

\EURUSDq J 

and here we can apply the martingale methods to yield 

Present Value = L'/^^^Notionabrf ( Q^^^^^^ - ' EURUSDt^ 



an expression that involves the second moment of spot. A model that has an infinite second moment 
for spot (and these things do appear in practice) will price such a contract at infinity. 

2. LIBOR paid in arrears. A similar type of deal occurs in the fixed income derivatives market 
under the name of LIBOR paid in arrears. Normally LIBOR is a rate fixed on a certain date Ti and 
paying at a later date T2. A simple contract depending on LIBOR is a ERA (Eorward Rate Agreement) 
whose payout is defined by 

S{L - K) paid at time T2 

here L is the LIBOR rate which fixes at Ti < T2 and 5 is the day-count fraction which is approximately 
T2 — Ti. If we use the T2-forward measure, the value of such a contract is simply the discounted 
expectation (5P(0, T2)'E{L-K) where P(0, T^) is the Ta-discount bond. Given that = | (P(t, Ti) - 
P{t, T2))/P{t, T2) is the ratio of a tradeable instrument by the numeraire, we know it is a martingale 
in the T^j -forward measure and since Lt^ is simply the LIBOR rate, the price of our ERA is 

P(0,T2)5E(L-ir) = P{Q,T2)5¥.{Lt,-K) (16) 
= P{^,T2)5{Lo-K) (17) 



where Lq = \ (P(0, Ti) - P(0, T2))/P(0, T2) is the so-called forward LIBOR rate, the rate that makes 
the ERA have zero value. A LIBOR in-arrears transaction is based on LIBOR paid at the wrong time, 
its payout being 

5{L - K) paid at time Ti 

A naive idea to price these transactions is to simply discount the by the forward discount factor 
P(0, T2) / P(0, Ti) which of course just affects equation (16) by a multiplicative factor and does not 
alter the fair forward price Lq. This is wrong because the expectation of L is no longer the forward 
LIBOR in the Ti-forward measure which this approach implicitly uses. Some banks have been arbi- 
traged in the past by using this naive approach. As in the case of the EURUSD note described above, 
if we convert the payment to a payment at time T2 then we can use the expectations in the T2-forward 
measure, the payout at time T2 is simply the accrued amount 



5(1 + 5L) -{L-K) paid at time T2 
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and its price will be 

P(0, T2)SE ((1 + SL) ■ {L - K)) 
and so the fair strike is strictly larger than the FRA rate 

Var(L) 

an expression that involves the second moment of the LIBOR rate. 

4.3 Dependence of the abscissae of convergence on the time. 

In this subsection we assume that = 0, so {St, t e [0,T]} is a martingale. In order to study the 
dependence of the abcissae of convergence on the time, we denote by u\ (t) and u*_ (t) the abcissae for 
Xf. The following proposition is a general property of a positive martingale. 

Proposition 4.1. Consider <t <t' <T. Then 

<(t)><(0 and u*_{t) <u*_{t'). 

Proof. 

Fix for a moment r > 1. The function (f){x) — on (0, oo) is convex. Since {St, t E [0, T]} is 
a positive martingale, assuming enough integrablity, t e [0, T]} is a submartingale. Hence, for 

t < t', 

E[Sl]<E[Sl,]. 
For e > 0, such that u*^{t') — £ > 1, we have that 

nsf^^'^-'] < ns;}^'^-'] < oo, 

and this implies u*^{t) > u*^{t') — e, and the result follows from (15). 

For the negative abscissa, just observe that for r < 0, the same function ip{x) — on (0, oo) is 
also convex, and apply the same reasoning. □ 

For Heston model, we can be more precise. From the bounds given in Theorem 3.4, it is deduced 
the behaviour of the abcissae for i ^ oo 



and for t ^ 0: 



lim u*_{t) — and lim u*_{t) — u^, 

t—*oo t—*oo 



limu*,(t) — oo and limu*_(t) — — oo. 



A plot of these funcions is given is Figure 3. 
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Figure 3. SoUd Une: Plot of u*^{t). Dashed line: Plot of u*_ (t). 



4.4 The effective vol-of-vol and the effective mean reversion 

It is interesting that in formula (10) the numbers u± only depend on the quotient u :— a/c. Expressed 
in terms of this parameter we have, for p 7^ ±1, 

l-2up± JAu'^ + 1 - Aujp 

= ^(T^) ■ 

Intuitively, this is a consequence of the fact that the parameter a (mean reversion) dampens the 
stochastic volatility whereas c (vol-of-vol) increases it, thus operating in opposite directions. However 
the parameters u± correspond to the moments at infinite time, and in general the relative strength of 
the parameters a and c in flattening the smile is time dependent. We propose to call this parameter the 
effective mean reversion factor, and its inverse c/a the effective vol-of-vol. The use of this parameter 
simpUfies the study of u± for p near ±1: 

lim U- = —- and lim u+ = 00, 

2u; + l 

and 

, if u; < 1/2 



— { ^ and \imu+ 

-00, ifu;>l/2 




5 Factorization of the moment generating function of Xf 

In this section we will work exclusively with the random variable Xt for t > fixed, and the time t will 
be considered a parameter. Denote by HC{a, 6, c, p, xq, vq, p, t) the law of Xt, that is, a probability 
on R that has the moment generating function given by (3). Observe that for all A > 0, 

7iC{Xa, Xb, Ac, p, xq, Xvq, p, t/X) = 7iC{a, b, c, p, xq, vq, p, t). 

In particular. 



HC{a, b, c, p, Xq, Vq, p, t) = HC{at/2, bt/2, ct/2, p, xq, vot/2, p, 2). 
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Since we are interested in a property true for all parameters, it suffices to prove that property for 
arbitrary a, 6, c, p, xq^vq, ji, and t — 2. So, in all proofs, we will take this value of t. 

In this section we use some powerful theorems of complex variable analysis that we apply to the 
complex moment generating function 



^{a-cpz)t/2 \ 2a6/c' 



cosh{P{z)t/2) + (a - cpz) s\n]i{P{z)t/2) / P{z) 



\_ - z') s\nh{P{z)t/2)/P{z) \ 

"""^^ \ cosh{P{z)t/2) + {a- cpz) smh{P{z)t/2)/P{z) J ' ^ ^ 

defined in a neighborhood of zero, introduced in Subsection 2.2. We will consider separately the 
different factors of this function, and at a later stage, we will combine them. 

5.1 The entire component 

Write 

F{z) = cosHP{z)t/2) + (a - cpz) ^^^h(P(^)t/2) ^^^^ 

Pyz) 

Recall that 

p{z) = [a — pcz)'^ + c^i^z — z'^) 

and P{z) = \/p{z). As in Subsection 3.2, but now in the complex plane, define the entire functions 
by the power series 



n=0 ^ ' 

and 

n=0 ^ ' 



We note that at each z. 



Lx{z) — cosh-\/z and Li{z) — smh^/z, 



independently of the branch of the square root. Indeed, in every neighborhood that does not include 
zero, the previous relations are true fixing an arbitrary branch of the square root. However, in the 
whole C, the functions Li and L2 are defined by the power series and not as a composition of cosh z 
or sinh z and a particular branch of ^/z. We consider the extension of F{z) to an entire function using 
the functions Li{z) and L2{z). Note that both Li{z) and L2{z) take real values on R, and that F{z) 
restricted to R coincides with the function F{u) defined in (13). 

The first interesting property of F{z) is that it has all the zeroes real and simple. This can be 
deduced from a deep theorem of Lucie [15, Theorem 1] who proves that all the singularities of the 
(complex) characteristic function of Xt are purely imaginary. In Appendix B there is an alternative 
proof of this property. Specifically, we prove that 



Proposition 5.1. The zeroes of F{z) are all real and simple. 
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Using the Hadamard representation theorem (see the Appendix, Theorem B.l), we deduce the 
following representation of F{z). 

Theorem 5.2. Let {an{t), n > 1} be the zeroes of F{z). Then 

oo 

F{z) = e-^/^e'^W^ 11 " ^Ti) { ^ I ' (^2) 

V an{t)J ^an{t)J 

where ^{t) e R. 

Remark 5.3. The number i/{t) in (22) is determinated by 



oo ^ 1 

F(l) = cosh (|a-pc|t/2)+sign{a-pc} sinh (|a-pc|t/2) = e'^^/^e'^^*) J]^ (^1 p-) expj— — }. 

(23) 



Remark 5.4. In all of this Section we exclude the case p = 1 and a = 2c because the moment 
generating function simplifies and has only one pole in that case. Then the factorization is trivial. 



5.2 The meromorphic component 

Now we will deal with the other factor of the function ^(z) given in (18), 

{z - z^)smh{P{z)t/2)/P{z) 



cosh{P{z)t/2) + {a- cpz) smh{P{z)t/2)/ P{z) ' 

Write 

_ (1 - ^) sinh(P(^)i/2)/P(^) 



cosh{P{z)t/2) + (a - cpz) smh{P{z)t/2)/P{z) ' 



where, as above, the function is extended to C using the entire functions Li{z) and L2{z). The 
function G{z) is meromorphic with simple poles at the roots of F{z). Using a theorem of Mittag- 
Leffler, see the Apendix, Theorem B.2, we prove 



Theorem 5.5. Le?{a„(t), n > 1} be the zeroes of F{z), ordered in the following way: < |ai(i)| < 
|a2(^)| < • • ■ • Then, for z ^ ttnit), \/n > 1, 

where bn{t) > 0, and^^ ^^j^ < oo, and the series converges uniformly in every compact set included 
in the disc \z\ < \ai\. 
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5.3 Factorization of the moment generating function of Xt 

We return to the real moment generating function M{u). Combining the results of the two preceding 
subsections we have that the moment generating function M{u) given in (3) can be factorized as (we 
abreviate a„(i) and &„(i) to a„ and 6„ for a moment), for w e (— |ai|, |ai|), 



M{u) = exp{du} 11(1 J exp <^ + vq 



2aVc' r 2ab u u%n 



n=l 

where d is a parameter to be defined later. Each exponential can be written as 

9nU 

CnU + -, . 

1 - u/an 

where 

-{vobn + 2ab/c^)/an and Qn^vobn/an- (25) 

Remark 5.6. When Aab/c^ = A; is a natural number, we can identify the moment generating function 
as the one given by Janson [13, Theorem 6.2], each eigenvalue with multiplicity k. That means, for 
such parameters, Xt is in the (non-homogeneous) second Wiener chaos. It is well known that for 
such combination of parameters the CIR model has the law of a sum of the squares of k independent 
Omstein-Uhlenbeck processes, and this fact has been used for many applications. See, for example, 
Grasselli and Hurd [10] and the references therein. 

On the other hand, for m in a neighborhood of 0, the moment generating function of each factor 
can be written as 

iV,M = e"(l-ii)-""^xp{-^}. 

\ On/ 11— U/an ^ 

where a„ • g'^ > 0. 

Proposition 5.7. For ^ > and C, 7 £ K, such that ( ■ ^ > 0, the function 



7VH = (l-V7)-'exp{^-^} 



/7- 

for u in a neighborhood ofO, is the moment generating function of an absolutely continuous random 
variable with probability density function 



(^^^ exp{-(x + C)7K2r(27A/C^)l(o,oo)(a;), ifC,l>0, 



or 



= exp{-(a; + C)7}^2r(-27\/Ca^)l(-oo,o)(a;), i/'C,7<0, 

where t — — l)/2 and is the Bessel function of index 2t. 
Proof. 

For 7 > 0, the moment generating function N{u) corresponds to the law of ^^1/(27) for a Bessel 
process 

Yt^C + '2^t + 2 f ^/YsdWs. 
Jo 
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See Revuz-Yor [17, Chapter 11]. Its density is also given in Revuz-Yor [17, page 441]. 

For < 0, consider the random variable Y defined above with parameters —( and —7, and let 
Y' — —Y. Its moment generating function is 



Its density is /iy (— x), where hy is the density of Y. And the result follows. □ 

Return to the factors Nj{u). The term e'^"" is a translation factor of the random variable considered 
above, so Nj (u) corresponds to a random variable with density hn{x — Cn), where hn is the probability 
density function given in Proposition 5.7. 

Finally, the product of characteristic functions of absolutely continuous laws with densities ki and 
/c2 corresponds to the convolution ki-kk^. 



Denote by i^'^^y^k.j the product ki-k - ■ - -kkn, that is defined without ambiguity because the convolution 
product is associative and commutative. 

With all these ingredients we construct a sequence of random variables that converges in law to Xt. 
In general, the convergence in law does not imply the convergence of the corresponding probability 
density functions. However, from a practical point of view, this fact does not matter. 

Theorem 5.8. The sequence of laws with densities -j^^^^hj converges to the law of Xt + d{t), where 

hn{x) = hn{x — On), and hn is the probability density function given in Proposition 5.7, with param- 
eters 





^ = 2a6/c^ 7 = an{t), C = 9n{t) 



where 



1. ttnif) are the roots of the function F{u). 



2. 



b„{t) 



4p(a„(t))(l -a„(t)) 



tp'{anit))p{an{t)) - 4:cpp{an{t)) -p'{an{t)){a - cpan{t)){{a - cpan{t))t + 2) 



3. Cnit) and Qnif) are given in equation (25). 



4. 



d{t) = Xq 



pabt 2ahv{t) vq{1 - e""*) 



c 2a 



5. ^{t) is given in formula (23). 
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Conclusions 

We have presented a new expression of the characteristic function of log-spot in Heston model that 
shows its good analytical properties and that facilitates its study. Through an analysis of the corre- 
sponding moment generating function, we give numerical formulas to obtain the abscissae of con- 
vergence, that have interesting applications. As examples we considered the computation of the 
parameters describing the asymptotic behaviour of the volatility smile for extreme strikes, and the 
verification that the model has enough moments to price wing dependent deals. Another application 
may be the possibility of assessing the stability of the moments in a calibration of a Heston model. 

In the second part of the paper, we factorized the moment generating function as an infinite product 
of Bessel type moment generating function. This gives a new insight of the Heston model, showing 
its complexity, and its relationship with Ornstein-Uhlenbeck processes. Further, each factor can be in- 
verted, and a sequence of random variables that converges in law to log-spot can be deduced. Though 
such sequence is not easy to manage, because it relies on the computation of the roots of a (real) func- 
tion and other parameters, and the convolution of densities involving Bessel functions. The fact that 
all computations are real can open the possibility of alternative methods to the numerical inversion of 
the characteristic function that are currently used by practitioners. 

Appendix 

Appendix A. Complex and real moment generating function 

For a sake of completeness, we recall some of the properties of the complex and real moment gener- 
ating function of an arbitrary random variable X. We follow Hoffmann- J0rgensen [12]. The complex 
function 

Mc(^)=E[e^^], 

defined on the set 

Dc^{zeC: E[|e^^|] < oo}, 

o 

is called the complex moment generating function (cmgf from now on) of X. We will denote by Dc 
the interior of D^. The restriction of Mc{z) to the real nimibers is the moment generating function 
(mgf) 

MM(u)=E[e"^], 

defined on the set 

L>R = {MeM: E[e"^] < oo}. 

o 

(we also put by the interior of D^). For the present purposes it is convenient to maintain the 
double notation with the subindices R and C. Since for 2; e C and a e R, | exp{a2;}| = exp{a ^z}, 
we deduce that 

Dc^{zeC: ^z e D^} ^ D^ + iR. 

The most important property of both Mc{z) and M^{u) is that they are analytic functions on the 
interior of its domains, and that the Taylor expansion of Mc{z) in a point of the real axis has the 
same (real) coefficients than M]^{u). Again, here, it is useful to introduce a double notation for the 
neighborhoods. For u e R and r > 0, we denote by -Br(m, r) the neighborhood centered at u with 
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radius r, and for z E C, Bc{z, r) is the neighborhood in C. If ^ {0}, the cmgf Mc (respectively 

o o o 

the mgf Mr) is analytic in L>c (resp. Dr). and for zq G-Dc there is r > such that 

^c(^) = \ (■^ - ■^o)"' ^ ^c(^o, r), (26) 

n=0 

o 

and for gDr there is r > such that 

Mm(«) = V ^ ^ i (k - iio)", Vii e SM(iio, r). (27) 

n=0 

o 

In particular, if G-Dr, then X has finite moments of all orders and, writing 

m„ = E{X% n > 1, 

we have that for some r > 0, 

oo 
n=l 

and 

oo 

Mr(k) = 5^^«", Vue5R(0,r). 

n=l 

Moreover, 

EfX'*] = M(")(0). 

Proposition A.l. Assume that there is a function $(^) of the complex variable z analytic in a neigh- 
borhood of such that 

for u in a (real) neighborhood ofO. Then the function ^{z) can be analytically continuated in a strip 
that includes the imaginary axis, and the characteristic function of X is 

Proof. Since M^{u) and Mc{z) have the same Taylor coefficients in a neighborhood of 0, we have 

o 

that ^{z) — Mc{z) in a complex neighborhood of 0. Since Mc{z) is analytic in Dc, the function $ 
can be continuated to that strip that includes the imaginary axis. □ 

Remember that is an interval of R (it may be 0), and the right extreme (respectively the left 
extreme) of Dr is called the right-abscissa (resp. the left-abscissa) of convergence. 

o 

Proposition A. 2. Assume that G-Dr, and that the right-abscissa of convergence, {3, (respectively 
the left-abscissa a) is finite. Then M-g^{u) has a singularity at (3 (resp. a). 

Proof. When X is positive, the cmgf coincides with the Laplace transform (except a change of sign 
on z) of the distribution function of X. By Widder [19, Theorem II.5b]), Mc{z) has a singularity at 
the real point (3. From the fact that in the real axis both Mc and Mr have the same coefficients of the 
Taylor expansion, Widder's proof can be translated to M^. For a general X, we can use the habitual 
technique of decomposing a bilateral Laplace transform as the addition of two unilateral ones (see 
Widder [19, page 237])). □ 
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Appendix B. Proofs 

Proof of Proposition 3.1 

We have that p(0) = and = [a — pY > 0, then, given the form or p{u), it follows that 
1 < u^. Hence, the affirmation 1] C D{Xt) is obtained from the other points of the proposition. 

1. Consider the case a > pc, p 7^ ±1. A bit of algebra shows that a — pcu- > and a — cpu+ > 0. 
So the straightline y = a — pcu is positive for u G m+] (see Figure 4). Since cosh x > 1, Va; e M, 
and sinha; > 0, Vx > 0, it follows that f{u) > in [«_,«_)_]. Note that if a = pc, then u+ = 1. 




Figure 4. Solid line: straightline y = a — pcu, case a > pc, for p > 0. Dashed line: 
parabola p{u) 

When p — —1, then U- = —a'^/{c{c + 2a)) < 0, and it is trivial that a + cu- > 0. So the 
straithgtline y = a + cuh positive for u > U-, and hence f{u) > in cxo). 

When p = 1, and c 7^ 2a, then = a^/ (c(2a — c)), and some calculations shows that in this case 
(a > c), then 

a a? 
> 



c ~ c(2a — c) ' 

and hence, a — cu+ > 0. Then a — cii > for all u < u+, and it follows f{u) > in (—00, it+]. 

2. Consider the case a < pc (note that this implies p > 0). Then a — cpu^ > and a — cp (^) =0. 
So the straighUne y — a — cpu is positive for u e a/(cp)]. Thus f{u) has no zeroes in such 
interval. On the other hand, the real zeroes of f{u) in u+] are the same of the ones corresponding 
to the function 

^^a-cp« 
P{u) 

because of the real hyperbolic cosine is never zero. For u > a/{pc), and the bound < tanhx < 1, 
for all X > 0, we have 

and when m e [a/ (pc), 1], from c^u{l — u) > it is deduced that 

a — cpu 
P[u) 

hence f{u) > also in [a/(pc), 1]. Now, we analyze the behaviour of f{u) for u e [1, u+\. Assume 
p ^ 1. From 

sinh (P{u)t/2) 

lim — 1. 

P{u)t/2 
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it follows that 



lim f{u) — 1 + {a — cpu+)t/2. 



Consider the straightline y = a — cpu. We have y{Q) = a and 7/(1) = a — cp < 0, and since u+ > 1, 
then a — cpu+ < 0. Let to = 2/(cp-u+ — a) > 0. 

(i) Ift< to, from {cpu+ — a)to/2 — 1, then for all u e [1, 

< {cpu - a)to/2 < 1, 



and thus 



Hence, Vm e [1, 



< {cpu - a)t/2 < 1. 



i^f^ - sinh (P(n)t/2) sinh (P(..)t/2) 

2 PHt/2 ^ PHt/2 ^cosn^nW^j, 

because smhx/x < coshx. Then f{u) > 0. 

(ii) Let t > to. Then /(I) > and /(«+) < 0, and thus there is at least one zero of f{u) in [1, «_,_]. 
The fact that there is only one zero is proved in Section 5. 

Finally, the case p = I and a < c or a = 2c are studied in a similar way. □ 



Proof of the results of Section 5. 

For easy reference, we recall here the two main theorems used in the proofs, expressed in the form 
that we need. The first one is the factorization Theorem of Hadamard (see, for example, Titchmarsch 
[18]) 

Theorem B.l. Let F{z) an entire function, P(0) ^ 0, with roots ai, a2, . . . , such that 

^ l/IOnT < 
n 

Then F can be represented as 

F{z) = P(0)e^" (l - — )e"/"". 

n 

The second theorem is due to Mittag-Leffler (Titchmarsch [18, page 110]) 

Theorem B. 2. Let G{z) be a meromorphic function such all the poles are simple, denoted by 
oi, a2, . . . , where < |ai| < \a2\ < ■ ■ • , and with residues at the poles bi, 62 • • • respectively. 
Suppose that there is a sequence of closed contours Cn such that Cn includes ai, . . . , a„ but no other 
poles, such that the minimum distance Rn ofCn to the origin tends to infinite with n, while the length 
ofCn is 0{Rn), and on Cn, f{z) — o{Rn). Then 

G(.).G(0)-.El^- 

n ' 
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In order to prove Proposition 5.1 we need three lemmas. 
Lemma B.3. For all 2; e C, 



coth V 2; I < sup I coth w I , 

\w\=J\^\ 



and similarly for tanh ^/z. 

Proof. 

Consider the formula 



, ,9 sinh^ a; + cos^ y 

coth^l^ = 2 2^- (28) 

sinh X + sin y 



Write z ~re^^ and choose an arbitrary branch of the square root, for example, take the principal one. 
Then 

I coth fz? - cos(^/2)) + cos2(v^ sin(^/2)) 

I V ^1 - g-^jj2^^ cos(^/2)) + sin2(yf sin(^/2)) 

sinh^(v^ cos (h) + cos^(v/r sin0) 

sup ^ ^ 

~ </>e[o,27r) sinh^(^ cos0) + sin^(v^ sin0) 

I 1 2 

= sup I coth w\ . □ 

\w\=^ 

Lemma B.4. For every e > there is no > 1 such that for n > no, 

sup I coth2;| < 1 + e 

|2|=(n+§)7r 

and 

sup I tanh2;| < 1 + e. 

\z\=mT 

Proof 

From the formula (28) it is clear that | coth z\ is the same for the points z — x + iy, —x + iy, x — 
iy, —x — iy. Hence, to bound coth z in the circle G C : \z\ = {n+ |)7r} we can restrict ourselves 
to study the arc with 3fj2; > i SJ^; > 0. (see the Figure 5 forn = 2). 



Figure 5. Arc \z\ = (n + |)7r, ?R:z>0,^z>0 for n = 2. 
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Given the periodicity coth(2; + kiri) — coth z, for all A; e Z, it suffices to bound the translation to 
the strip < Qz < n (see Figure 6 (a)). For K big enough, for n > uq all translations are included 
in the shaded region Ak U Bk of Figure 6 (b). Thus, it suffices to prove that for K big enough, 

^^PzeAkuBk I z\ < 1 + e. For z = x + iy e Ak, (that is, x > K), we have 



I coth ^ I 



sinh^ X + cos^ y ^ sinh^ x 



sinh^ X + sin^ y sinh^ x 
that goes to 1 when K ^ oo. 



1 + 



sinh X 



2 „ — 



< 1 + 



sinh^X' 





(a) 



(b) 



Figure 6. (a) Translation to the arc of Figure 5 to the strip < < tt. (b) The 
translation for all n > no are included in the shadow region 

Now we study the bound on the triangle Bk of the Figure 6 (b). By the principle of the maximum, 
we need only to study the function on the border of Bk- On the vertical side, z = K+iy, y e [0, 7r/2], 
works the bound that we have found above. On the horizontal side, z = x + i^, x e [0, /T], by the 
formula (28), 

siiilr X 



coth 2; I 



< 1. 



sinh^ X + 1 

Finally, the hypotenuse is z ~ x + iy with y e [0, 7r/2] ix — —2Ky/7r + K. For y e [7r/4, 7r/2], we 
have that cos 2y < 0, and by 

2 cosh 2a; + cos 2|/ 



coth 2;! 



cosh 2x — cos 2y ' 



we deduce that 



I coth^l < 1. 

For y e [0, 7r/4], we have x G [K/2, K], and again by (28) we obtain the bound. 

To bound | tanh 2; | for — nn we do the same reductions as before, and use the relationship 



tanh z = coth (2; — i— ) , 

to translate the bounds of coth z to tanh z. □ 

We also need the properties of the contour determinated by the polynomial j^iz). For d > let 
be the contour 

and denote by its length and by R^, the minimum distance from to the origin. See Figure 7 
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Lemma B.5. With the above notations, for d big enough, is a homotopic to 0, lim^^oo Rd — 
and Ld ~ o{R'^), when d — > oo. 

Proof 

Since p{z) is a second degree polynomial, there are two constants Ki, K2 > such that for \z\ 
big enough, 

Ki\z\^ < \p{z)\ < K2\z\^. 

It follows that for d big enough, the circles with radius d/K2 i d/Ki bound lower and upper the 
contour Cd- This implies that Cd is closed. 




Figure 7. Contour \p{z)\ — d, for d big. 



The polynomial p{z) has real coefficients and has a positive and a negative root (p 7^ ibl), and 
since the length of a contour does not change by translation, we can consider that p{z) = z'^ — a, with 
a > (change also d in order that the coefficient of z^ is 1). Write z — x + iy; then the equation that 
determines C„ is 

{x^ + y^ + aY - Aax^ = 
In polar coordinates, x — r{9) cos 9 and x — r{9) sin 9, the contour is given by 

r^{9) - 2a cos{29y{9) + a^-d'^, 9e [0, 2tt]. 

For d big enough, we need only to consider the solution 

r{9) = \Jacos{29) + y/a^cos{A9) 
that determines a curve clearly homotopic to 0. To compute the length Ld, by symmetry, 

= 4 Jr\9) + {r'{9)Yd9. 
Jo 

We have limd-^oof{d)/d — and liuid-^oo f' / d — 0, uniformly in ^ e [0, 7r/2], and the lemma 
follows. □ 

Proof of Proposition 5.1 

Remember that we consider t — 2. First, note that for p 7^ ±1, 

lim - 'P'^" = 

z^oo (a — cpzy + c^lz — z"^) p2 — 1 
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Hence, 



lim 

2— »0O 



a — cpz 



P{z) 



(29) 



Consider the following three cases: 
Case 1. \p\/^/l-p^ > 1. 

First step. The objective of this step is to prove that in the contour 



for n big, 



Cn^{zeC: \p{z)\ 
coshPf^;)! < {a — cpz) 



n + 0.5)V}, 
sinh P(z) 



P{z) 



(30) 



Then, by Rouche Theorem, both functions F{z) and (a — cpz) ?,m\i{P{z)) / P{z) have the same num- 
ber of zeroes in C„. To prove the inequality (30), take e > Q such that I < 1 + e < \p\/ \/l — p^. By 
Lemmas B3 and B4, for n big enough, and z G C„, we have | coth P{z) \ < 1 + £, and apply the limit 
(29). 

Second step. Study of the zeroes of F{z). Assume that p 7^ 0. The zeroes of [a— cpz) smk{P{z))/ P{z) 
are z — a/{cp) and the zeroes of smh{P{z))/ P{z). We have that 



sinhP(z) 



n=l 



P{z) 



Hence, the zeroes of this function are the roots of the second order polynomial 

p{z) — — n^TT^, n — 1, 2, . . . , 

that we denote by a±n; write also u- — a-o and u+ — a+o, in agreement with the notations of 
Section 3. They all are real and 

• • • Q;-(n+i) < tt-n < • • • a-1 < a-0 < < a+o < a+i < • • • < «+„ < a+(n+i) < • • • ■ 

So, for n big, F{z) has 2n + l zeroes in C„. 

Now we count the number of real roots of F(z) in C„. Assume that a > cp and p > 0. We saw 
in the proof of Proposition 3.1 that a/ (cp) > and assume that a/ (cp) e (a+fc, Q;+(fc+i)) (the case 
a/{cp) = a+k needs to be studied as a particular case). Then, 

(i) In each interval («+„, Q;+(n+i)), forn = 0, . . . /c — 1, F{u) = f{u), where f{u) was defined in 
(12), and /(u) has one root. This is deduced because the roots of f{u) in such intervals are the 
solutions of _ 

P{u) 



iaj\P{u) — 



a — cpu 



(ii) In (a+fc, a+(fe+i)), f{u) has 2 roots. This claim is proved, observing that a — cpa+k > and 
a — cp a+k < 0, and the the curve —P{u) /{a — cpu) cuts two times the curve tan P{u) in that 
interval. 
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(iii) In each («_(„+!), a_„), n > 0, the function f{u) has one root. This is proved as in point (i). 

All the other possibilities for a, c and p are discussed in a similar way, and we obtain that F{z) has at 
least 2n + 1 real roots in C„. So the Theorem follows. 



Case 2. \p\/ y^l — < 1. Here use the contours 

C; = {^eC: b(^)|=nV}, 

and prove that for z e C'^,(n big), 

sinh P{z) 



[a — cpz)- 



P{z) 



< coshPO^ . (31) 



and finish the proof as in Case 1. 



Case 3. p — ±\/2/2. Write p„ = p+l/n, for n > 4, and let F„ be the function F with p changed 

by Pn- We have Fk{z) F{z) as A; ^ oo, uniformly in every disc. By Hurwitz theorem (see 
Titchmarsch [18]), the roots of F{z) in such disc are the limit points of the roots of Fk{z) in the disc. 
So the roots of F{z) are also real. □. 

Proof of Theorem 5.2 

In the proof of Proposition 5.1 we saw that for n big, positive or negative, in each interval 
(«+„, q;+(„+i)) and (q;_(„+i), q;_„) there is one and only one root of F[z), where 

p{oi±n) = -nV. 

Since p{u) is a second degree polynomial, there are constants Ki, K2 > Q such that for \u\ big 
enough, 

.2 ^ \^f„.\ \ ^ 1^ „.2 



and then 



Kxu' < \p{u)\ <K2u\ 



nn < \a±n\ < nir. (32) 



Hence, 



n "■ n ' ' 

By Hadamard factorization Theorem B.l we obtain the representation (22). □ 

Remark B.6. An alternative way to prove the previous result is using that the order (as entire function) 
of F{z) is less or equal to 1. This can be deduced from that fact that by the definitions (20) and (21) 



|L,(^)|<cosh(Vkl), J = l,2, 
and that the order of Li{z) can be easily computed from (20) and is 1/2. 
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Proof of Theorem 5.5 

In a similar way that in the proof of Proposition 5. 1, we are going to prove that there is a constant 
C > such that for z E Cn or z E C!^, for n big enough, 

\G{z)\<C, 

where C„ and C'^ where defined in Theorem 5.1. Then, thanks to Lemma B.5, we can apply the 
theorem of Mittag-Leffler B. 2 that gives the expression (24). Consider three cases: 

Case 1. |p|/\/l - p2 > 1. Take e> such that 1 < 1 + £ < \p\/\/l - p^. By Lemmas B3 and B4, 
for n big enough, and z e Cn, \ coth P{z) | < 1 + £. Therefore 



\G{z)\ 



(1 - .:)/P(.:) 



< 



cothP(^) + (a - cpz)/P{z) 

\{l-z)lP{z)\ 
\{a-cpz)/P{z) \ - \ coi\iP{z)\ 





(1 








cothP( 






{a-cpz)/P{z)\ 



< c. 



Case 2. \p\/\/l - p^ < 1. Let 5 > be such that \p\/\/l - < 5 < I, and £ > such that 
(1 + e)5 < 1. Then, for z E Cnin large). 



(1-^) tanhP(^)/P(^) 


< 


tanhP(^) (1 - 


z)/P{z)\ 


1 + (a - cpz) tanh P{z) /P{z) 


- 1 


- tanhP(z) (a- 


-cpz)/P{z)\ 



In both cases 1 and 2, by Mittag-Leffler Theorem B. 2, 



°° h 



n=l 



where is the residue of G{z) in the pole a„. Since G{z) is the quotient of two entire functions, and 
the pole is simple, 

(1 - an) sinh(P(a„))/P(a„) 



F'iar. 



using that aj is a root of F{z), differentiating and simpUfying we obtain that the corresponding residue 
is 

^ ^ 2j;(a,,)(l - (i„) 

p'{an)p{an) - 2cpp{an) - p'{an){a - cpan){a - cpun + 1) ' 
Consider the function defined by the previous relation: 

( ) — 2p(u)(l — u) ^^^^ 

p'{u)p{u) — 2cpp{u) — p'{u){a — cpu){a — cpu + 1) 



Then 

ifpVl 

u-»±oo"^^ ' if = 1 and c 7^ 2ap. 



lirn ^(i.) = <( f ./^ _ (35) 
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So, for large n, it is clear that fo„ > 0. For small n the positivity of is proved through an analysis 
of the sign of f{u) and f'{u) in the different intervals where there are located the roots of /. These 
roots and its location are clearly studied in Lucie [15]. 

Case 3. p — ±V2/2. As in case 3 of Theorem 5.2, the result is obtained by continuity, using 
monotone convergence Theorem. □ 
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